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SUMMARY 


Experimental results are presented of an Investigation on the flutter 
of propeller blades with three different plan forms. Data for a proposed 
supersonic propeller having highly twisted 2-percent -thick blades Indicate 
that the flutter boundaries are complex and that several modes of flutter 
are Involved. Data for a 6-percent-thick blade far which, both the high 
Mach number and the low Mach number flutter boundaries are obtained In 
detail Indicate that Mach number effects can be beneficial In alleviating 
flutter. Results for two 3-percent-thick blades having different Internal 
damping show that an Increase In the Internal, damping 1 b beneficial in 
reducing the magnitude of the stresses both during flutter and for condi- 
tions of blade resonance. 


INTRODUCTION 


Stall flutter Is recognized as one of the most serious vibration 
problems encountered In propeller and compressor operation. Although a 
large amount of information Is available from full-scale tests over a 
period of time from World War I to the present, these data are difficult 
to corr elat e because of the wide variations of the test conditions. A 
few systematic studies are available In which attempts were made to Iso- 
late the effects of some of the significant parameters In the problem. 
Reference 1, In particular, contains data on the effects of torsional 
stiffness, section thickness ratio, sweepback, length-chord ratio, section 
center-of -gravity location, blade twist, blade taper, Mach number, and 
fluid density. A criterion is also presented In reference 1 which stipu- 
lates the conditions necessary for flutter-free operation of propeller 
blades. 
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The present paper contains results of an Investigation an the flut- 
ter of propeller blades with three different plan farms. Some of the 
data were obtained In the same wArmey as those of reference 1 and are 
extensions of that work. Specific attention was given to tests of a sug- 
gested supersonic propeller having highly twisted 2-per cent -thick blades. 
This model was operated over a typical blade-angle range in order to 
def ine its boundaries far torsion, bending, and bending-torsion modes of 
flutter. lests were also made with a 6-percent -thick propeller In which 
the speed of sound of the medium was changed In order to study the shape 
of the flutter boundaries, particularly for the higher Mach numbers of 
operation. In addition, some effects of structural damping on the flut- 
ter boundaries as well as on the blade stressed during operation near the 
flutter boundary were studied for an untwisted 3-percerrt-thick propeller. 


SYMBOLS 


a 

b 

D 

El 

f 


f hl 




«h 

GJ 


speed of sound In operating medium, ft /sec 
blade semichord, ft 
sound-speed coefficient 

blade diameter, ft 

bending stiffness, lb-ft 2 

frequency of recorded oscillation, cps 

natural first bending frequency, cps 

natural second bending frequency, cps 

natural frequency at propeller rotational speed, cps 

natural frequency at zero rotational speed, cps 

natural torsional frequency, cps 

torsional damping coefficient 

bending damping coefficient 

torsional stiffness, lb-ft 2 
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V 

V /^ 

P 

V 
P 

“a 


blade thickness, ft 

polar moment of Inertia about section center-of -gravity- 
location per unit length, slug-ft 2 /ft 

Southwell coefficient 

blade length, ft 

Mach number 

blade mass per unit length, slugs /ft 
propeller rotational speed, rps 

non d Inte nsions! radius of gyration of blade section about sec- 
tion center-of -gravity location, ]j I C g/mb 2 
blade velocity, ft /sec 
flutter-speed coefficient 
blade angle, deg 
relative-density factor 

mass density of operating medium, slugs /cu ft 
natural torsional circular frequency, radians /sec 


Subscript : 

0.81 referred to 0.8-blade -length station 


APPARATUS 
Test Facility 


Tests were conducted in the Langley vac uum sphere, previously 
described in reference 1, in which both the speed of sound and the pres- 
sure of the medium were varied. Pressures ranging from l/6 atmosphere 
to full atmosphere were used in the present tests. Mixtures of Freon-12 
and air were used in various proportions to give speeds of sound from 
approximately 550 to 1,100 feet per second. A 500-horsepower variable- 
speed electric motor was used to drive the propeller Tnndpi h which were 
tested at zero advance except far the Induced flow. 
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Model Blades and Hubs 

The plan forms of the four model blades tested are shown In figure 1 
and additional Information about the physical characteristics of each 
blade Is given In table I. These models were all tested as one-blade 
configurations with suitable counterweights attached as in reference 1. 

Model 1 (model of proposed steel supersonic propeller) was a twisted 
blade with NACA 16 -series symmetrical airfoil sections and was of solid 
aluminum-alloy construction. It was tested in a fixed-pitch hub for the 
purpose of defining the flutter regions for an extremely thin ( 2-percent - 
thick) blade. Blade -plan-form curves for this model are given in figure 2. 

Model 2 had an NACA l6-006 airfoil section, was untwisted, and was of 
solid aluminum— alloy construction. Thl « blade was tested In both the 
fixed- and controllable-pitch hubs for comparison. Most of the testing 
was done with the controllable-pitch hub for the purpose of mare carefully 
determin in g the boundaries of the flutter regions, particularly the upper 
boundaries. 

Models 3b. and 3b were Identical. In plan, farm and differed only In the 
method of construction. They were untwisted and had NACA 16-003 airfoil 
sections. Model 3a was solid steel whereas model 3b consisted of two 
laminations of steel bonded together flatwise with Cycleweld cement. 

The objective In laminating model 3b was to Introduce additional Internal 
blade damping In the same manner as described in reference 2. Both of 
these models were tested In a fixed-pitch hub. 


Instrumentation 

Bending and torsion strains near the root of the blade were obtained 
by means of strain gages located as shown in figure 1. Oscillograph 
recordings were made from which frequencies of oscillation were deter- 
mined. The rotational speed, which for zero forward velocity is equiva- 
lent to the resultant velocity, was recorded on the same record. Blade 
angles were measured while the propeller was at rest in all cases. Twist 
angles were measured optically with the aid' of stroboscopic light during 
the tests of model 1 and have been Included In the presentation of the 
data as an aid In the Interpretation of the test results. 


TEST PROCEDURES 


The same general* procedures were used In the present tests as for 
the tests of reference 1. At a given atmospheric pressure and blade 
angle, runs were made at various rotational speeds until flutter was 
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obtained. Results of several runs at different blade angles were useful 
in defining the flutter boundaries. For seme cases In which a controllable- 
pitch hub was used, both blade angle and rotational speed were varied during 
the run to expedite the definition, of certain portions of the flutter 
boundary. 

Flutter speed, as defined In this paper, is the lowest speed at which 
a sustained oscillation in either bending or torsion was detected. In 
some Instances for models 3 & and 3 b where flutter was not severe, the flut- 
ter speed was exceeded in order to obtain stress data inside the flutter 
boundary. It was also possible. In same cases, to pass through a given 
flutter region and to define other flutter regions at higher speeds. 

Damping values as listed In the paper were determined from the decre- 
ments of strain-gage signals recorded while the blades were mounted In the 
hubs, which were In all cases mounted on the motor shaft, but not rotating. 


RESORTS AND DISCUSSION 


The results of flutter tests on the propeller models listed In 
table I are discussed under two headings. Hie effects of such parameters 
as thickness, built-in twist, Mach number, and structural damping on the 
shape of the flutter boundaries are first discussed. These data are pre- 
sented conventionally In the form of flutter-speed coefficients as a func- 
tion of blade angle. Same discussion is then given to the bending and 
torsional stresses measured in the vicinity of the stall flutter boundary 
for two models differing in their amounts of Internal blade damping. 


Effects of Various Parameters an the Stall 
Flutter Boundaries 

Thif-knftas and built-in twist .- The results of the flutter tests on 
model 1 are presented In figure 3 . Data are plotted in the farm of 
flutter-speed coefficients as a function of 'blade angle (which, as 
explained later. Includes blade twist) for density values of the operating 
medium of 0.00079 and 0.00040 slug/cu ft. In the term , V 

Is the velocity of the propeller blade at the 0 . 8 -blade -length station, 
b Is the semichord at that station, and Is the blade natural tor- 
sional circular frequency. Also noted In the figure are the flutter 
frequencies which can be compared with the rotational frequencies of the 
right-hand vertical scale. In general, the recorded frequencies of the 
data plotted In figure 3 cure not Integral multiples of the rotational 
speed. 
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Model 1 was tested at reduced pressure to make the present test con- 
ditions comparable to the operating conditions of the proposed full-scale 
propeller. During these tests, the model was operated -over a wide blade- 
angle range In order to define Its flutter boundary. Because this blade 
was only 2 percent thick and because of a fairly large amount of built-in 
twist, particular attention was given to a determination of the amount of 
twist-up encountered during testing. The twist angles were measured 
optically for the 0 . 8-blade -length station and have been Indicated by the 
lines of long dashes In figure 3 as an aid In the Interpretation of the 
data. The forces on the blade tending to twist it result from the fol- 
lowing components: (l) the aerodynamic pitching moments, (2) the chard- 

wise components of the centrifugal forces (ref. 3) , and (3) the radial 
components of the centrifugal, forces tending to straighten the blade 
(ref. it)* AH these effects are additive In the negative blade-angle 
range and, hence, the amount of twist -up is large. For positive blade 
angles, components (2) and ( 3 ) oppose each other and, hence, the amount 
of twist -up is less. 

The tests were conducted by Initially setting a given blade angle 
and then Increasing the speed until flutter was encountered. The long- 
dashed lines In figure 3 are In a sense operating lines since they Indi- 
cate the actual blade angle aa measured during the run at various speeds 
as opposed to the static blade angle measured at the start of the run. 

In some cases. It was possible to pass completely through one or more 
flutter regions before the sustained oscillations became sufficiently 
violent to curtail the testing. As Indicated In the figure, three types 
of flutter were obtained for this model. In general, torsional oscilla- 
tions and coupled bending-torsion oscillations were associated with 
milder flutter than the bending oscillations. As a matter of Interest, 
the blade failed while operating near the bending flutter boundary. 

The data of figure 3(a) Indicate that it would, be very difficult to 
schedule the operation of this propeller In such a manner as to avoid 
these flutter regions. At a lower ambient pressure, however, as In fig- 
ure 3(b ) , the flutter boundaries are somewhat different and It appears 
that the propeller could, be Increased In speed without encountering flut- 
ter, provided the blade angle were changed with speed to avoid the flutter 
regions. The bending flutter region was eliminated entirely in the range 
of these tests, and the bending-torsion flutter region was greatly reduced 
In extent. The torsional oscillations encountered were also mild for the 
range of blade angles tested and were safely traversed. No attempt was 
made at this value of the density to evaluate the minimum flutter-speed 
coefficient which would normally be associated with the higher positive 
blade angles. 

Mach number .- There was some indication In the work of reference 1 
that the flutter regions of propellers might be limited mainly to the 
subsonic tip-speed range. Ulus, at a given blade angle, if flutter were 
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not encountered up to a tip speed corresponding to a Mach, number of 
approximately 1.0, no flutter would be encountered at M ghft-n tip speeds. 
Shis result suggested that the upper extent of the flutter region was 
limited because of Mach number effects. Shis phftwrHnpnrm is Illust rated 
schematically In the following sketch where a hypothetical flutter region 
having both an upper and a lower boundary is plotted: 



If operation of the propeller at point 0 were desired, a schedule of 
operation In which the blade angle was first Increased to A and than 
the tip speed was Increased to the desired operating point might result 
In destruction of the blade since part of the flutter region would be 
traversed. On the other hand if the tip speed were first Increased to 
B, the blade a ng l e could then be safely Increased to the desired operating 
point without encountering flutter. She use of a controllable-pitch hnh 
In these tests of model 2 had the objectives of ni ng both the upper 
and lower flutter boundaries and of exploring in more detail the possi- 
bility of this type of operating procedure. 


She main variable In these tests was the speed of sound of the 
operating medium. Results are shown In figure 4 In which flutter-speed 
coefficients are plotted as a function of blade an gi#» far three values 
of the sound-speed coefficient. Sound-speed coefficient is dafin^d as 

(bcD^/a) , where b Is the blade semichord at the 0.8-blade-length 
0 .0 "L 


station, m a Is the first torsional circular frequency of the blade, 

and a Is the speed of sound of the operating medium. At the lowest 
value of the sound-speed coefficient, which 1 nc. i fl eirhai i y corresponds to 
the full atmospheric density In air for this model, only the lower por- 
tion of the flutter boundary was defined because of stress limitations. 

In general, it can be seen that the flutter region becomes smaller as 
the sound-speed coefficient becomes larger. The trend of figure 4 indi- 
cates that there would be no flutter region at an at a value slightly 
hi gher than 0.31. This trend Is a confirmation of the findings of refer- 
ence 1 and appears to be a Mach number effect. 
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The data of figure 4 are replotted In figure 3 to show this effect 
more directly as a function of Mach number. The Mach number at flutter 
Is plotted as a function of blade angle far three values of the sound- 
speed coefficient. It can be seen that flutter occurs over a «Tna.nm» 
blade-angle range at the higher Mach numbers . At a sufficiently high 
Mach number, safe operation of this propeller above the flutter boundary 
Is possible. The Mach number was first Increased at low blade angles 
and then the blade angle was Increased to higher values while Mach num- 
ber remained constant. In a few cases, it was then possible to reduce 
Mach number for the purpose of defining the upper flutter boundary. Far 
this particular propeller, the classical-flutter speed was either higher 
than could be obtained In the tests or did not exist at all far these' 
conditions . 

Structural damping .- Increasing the structural damping has been 
suggested as a way of Increasing the minimum stall flutter speed of pro- 
peller blades. One method that suggests itself is to add damping at the 
root of the blade. Another method is the addition of damping In the 
blade itself as has been suggested In reference 2. Some attempt has been 
made to vary both the blade and hub damping during the present teats and 
the results of these tests are discussed In this section. 

The effect of hub damping on stall flutter is shown In figure 6 
where the flutter-speed coefficient Is plotted as a function of blade 
angle for two values of the torsional damping coefficient. These data 
are for model 2 tested in full atmospheric density. The structural 
damping coefficients apply to the hub and blade In combination and were 
obtained -while the blade was not rotating. The difference In the damping 
la believed to be due to the differences In the hubs since the same blade 
la used in both cases. The results of figure 6 indicate that Increasing 
the torsional damping coefficient from g a = 0.004- to 8^ = 0.022 was 

effective In increasing the minimum flutter-speed coefficient by approxi- 
mately to percent. Ibis result la in conformity with those obtained far 
the thin wings In reference 5- 

In order to study the effect of Internal blade damping, models 3a 
and 3b were fabricated and tested. Model 3a was of solid steel construc- 
tion and thus had relatively low Internal damping. Model 3b was of 
laminated construction designed to give increased Internal damping. The 
method of construction was the same as for model V of reference 2 inas- 
much as that model seemed to give the largest Increase in Internal damping 
over the basic configuration. 

The values of torsional damping coefficient for each of these 

two models mounted In a fixed-pitch hub on the motor shaft are given In 
figure 7 as a function of vibratory torsional stress. The damping in 
each case 1 b seen to be approximately a linear function of the stress 
far the range of the tests. For the solid blade, the damping Is seen to 
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vary from approximately 0.0023 to 0.00*4-3* "which is of* the same order of 
magnitude as the values obtained In reference 2 during .bench tests. The 
overall damping values are seen to be approximately 50 percent higher 
far the laminated blade than for the solid blade at comparable stress 
values. It is believed that this Increase over the basic configuration 
is an indication of the change in the blade damping since both models 
vere tested in the same hub arrangement. This Increase in damping is 
much less than had been anticipated on the basis of results in refer- 
ence 2. 

She effects an the flutter boundaries of the additional damping 
obtained with model 3b are shown in figure 8. She flutter-speed coeffi- 
cients are plotted as a function of blade angle for two different densi- 
ties. At full atmospheric density (p = 0.00238 slug/cu ft), there are 
no noticeable differences in the flutter boundaries far the range of the 
tests; however, at l/k atmospheric density (p = 0.00039 slug/cu ft), the 
blade having the higher damping also had somewhat higher flutter-speed 
coefficients as Indicated in figure 8(b). 


Blade Stresses Bear the Flutter Boundary 

Some stress data were recorded far models 3a and 3b in order to 
evaluate the effects of the additional interna], damping of model 3b. It 
was noted that in the vicinity of the flutter boundary several resonance 
conditions were encountered during which rather high stresses were noted. 
In order to assist in the interpretation of some of the stress data, 
figure 9 has been prepared to give the frequency spectra of these two 
model blades. 

The frequency spectra shown in figure 9 by the solid lines consist 
of the first bending, second bending, and first torsional blade fre- 
quencies plotted as a function of the rotational speed. Also shown by 
the dashed lines are the rotational, harmonics of order 1 to 9* The 
values of the blade frequencies at zero rotational speed are measured 
and are approximately the same for both models. The effects of rotational 
speed on the first and second bending frequencies were calculated by means 
of the Southwell relation: 


(M 2 ‘ ( f n»o) 2 + El 2 

where K is the Southwell coefficient which depends fpr its value an the 
blade geometry, angle of attack, and the mode of vibration. For the cal- 
culations of this paper, the values of K far the bending frequencies 
were deter mine d by extrapolating those presented in reference 6 to the 
blade-offset values of these tests. The effect of rotational speed on the 
torsional frequency was estimated from the recorded oscillation frequencies 
for a blade angle of 2*4-°. 
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The blade torsional and bending vibratory stresses, respectively, 
are shown In figures 10 and 11 as a function of rotational, speed far 
models 3 s - and 3b > Data were obtained for two blade-angle settings near 
the stall and for air densities of 0.00238 and 0.00039 slug/cu ft. Data 
shown in the figures are considered representative of those of the tests 
and allow a comparison of the stresses directly for two blades which have 
different Internal damping but which are otherwise nearly Identical. It 
is also possible from the data presented to note the effects of density 
of the medium on the measured stresses. 

During the tests, several major resonance points were noted for 
which the torsional, stresses were fairly high. The maximum stresses at 
the rotational, speeds of these resonances were recorded and are presented 
In figure 10. The stresses are gene-rally lower for the blade that has 
the larger amount of damping (model 3b). As was noted previously the 
sharp rise In the stresses at flutter was not appreciably retarded at 
atmospheric density of 0.00238 slug/cu ft as Indicated in figures 10 (a) 
and 10(b). At a lower atmospheric density (p = 0.00039 slug/cu ft), the 
rise In flutter stresses is delayed somewhat for the blade with addi- 
tional damping. At the higher density, it was not possible to penetrate 
the flutter region to any great extent by Increasing speed because of 
this rapid rise In the stresses. At the lower density, however. It was 
possible to penetrate a considerable distance Into the flutter region 
before limiting stresses were reached. As a matter of Interest, for the 
conditions of figure 10 (d) , It Is believed that no flutter was obtained 
although rather high vibratory stresses were recorded at seme of the 
resonances encountered. Although decreasing density has a beneficial 
effect on the flutter stresses, those stresses measured during resonances 
axe still, severe and of the same order of magnitude as the flutter 
stresses. 

He se latter data, which Illustrate the response of a low-damped sys- 
tem to a periodic forcing function, point up the problem of Interpretation 
of data obtained In such tests. In the case of resonance peaks such as 
are shown In figure 10(d), the frequency of oscillations tend to he Inte- 
gral multiples of the rotational speed; whereas. In a flutter condition, 
this Is not normally the case. It Is thus convenient to have a spectrum 
plot of the type shown in figure 9 as an aid In Interpretation of results. 

The vibratory bending stresses, which were obtained at the same time 
as the vibratory torsional stresses just discussed, are shewn in figure 11. 
Several major resonance points occur also In the bending records. These 
major responses are found to occur when the natural frequencies of the 
blades coincide approximately with Integral multiples of the propeller 
rotational speed. In general, the bending stresses tend to he lower than 
the torsional stresses In this operating range. Also, the flutter mode 
is primarily a torsional frequency In this operating range since the 
bending stresses do not clearly define the condition for the onset of 
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flutter. As in the case of the torsional stresses, the "bending stresses 
generally show a decrease with decreasing density far the . speed range of 
the tests. At the lower density, there is also a tendency for the meas- 
ured "bending stresses to "be lower for the blade with the higher damping. 


COHCHJBIONS 


An investigation relating to the flutter of thin propeller blades 
with three different plan forms has been conducted and the following 
results have been obtained: 

1. A highly twisted 2-percent -thick blade was observed to have 
rather complex flutter boundaries and to flutter in bending, torsion, 
or bending-torsion modes depending on its operating conditions. 

2. Both the upper and lower flutter boundaries were defined far a 
6-percent -thick blade for which the Mach number effects appeared to limit 
the extent of the flutter region. This phenomenon made possible flutter- 
free operation at tip speeds corresponding to a Mach number of approxi- 
mately 1.0 far all blade angles. 

3. Increases in the mi.nl.rmnn flutter-speed coefficients were noted 
for increases in structural damping and, in particular, hub damping. 

4. Internal blade damping was shown to be generally beneficial in 
reducing the magnitude of the vibratory stresses of a 3-percent-thick 
blade during flutter. 

3. At some resonance conditions, stresses of the same order of mag- 
nitude as the flutter stresses were encountered. These stresses were 
also lower for Increased internal blade damping. 


Langley Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , March 28, 1956. 
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(a) Model 1 (model of proposed supersonic propeller). Hub radius, 
1.25 feet; blade thickness, 2 percent; twisted 31 ° from root to 
tip (see fig. 2). 



(b) Model 2 used for Mach number studies. Hub radius, 1.27 feet 
blade thickness, 6 percent; zero twist. 



(c) Models 3a and 3b used for structural damping and stress studies. 
Hub radius, 1.225 feet; blade thickness, 3 percent; zero twist. 

Figure 1.- Propeller flutter-test blades. 
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Blade length , percent 

Figure 2.- Blade-plan-form curves far model 1. 
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(b) p « 0.00040 slug/cu ft 
Figure 3 .- Concluded. 
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Figure If.- Flutter-speed coefficient as a function of blade angle 
for three values of the sound-speed coefficient. Model 2 
p *■ 0.00258 slug/cu ft. 
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Figure 5-- Mach number at flutter as a function of blade 
far three values of the sound-speed coefficient . Model 2 
p ■» 0.00238 slug/cu ft. 








Figure 7.- Variations of 
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% Figure 9- - Frequency spectra for models 3a and 3b. 
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